Exact ground states of quantum spin-2 models on the hexagonal lattice 
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We construct exact non-trivial ground states of spin-2 quantum antiferromagnets on the hexagonal 
lattice. Using the optimum ground state approach we determine the ground state in different 
subspaces of a general spin-2 Hamiltonian consistent with some realistic symmetries. These states, 
which are not of simple product form, depend on two free parameters and can be shown to be 
only weakly degenerate. We find ground states with different types of magnetic order, i.e. a weak 
antiferromagnet with finite sublattice magnetization and a weak ferromagnet with ferrimagnetic 
order. For the latter it is argued that a quantum phase transition occurs within the solvable 
subspace. 

PACS numbers: 75.10.Jm, 73.43.Nq, 75.25.-fz 
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I. INTRODUCTION 

Usually one faces the problem of finding the ground 
state and excitations of a given (many-body) Hamilto- 
nian that is supposed to describe the physical problem of 
interest. The solution then allows, at least in principle, 
to determine other quantities, e.g. correlation functions. 
Usually an exact solution is not possible, especially in 
dimensions d > 2, and one has to rely on approximations 
(e.g. perturbation theory) or numerical methods. 

Another approach does just the opposite 1 . It starts 
with a given ground state \ipo) and tries to find the cor- 
responding Hamiltonian. This can be achieved by find- 
ing positive-semidefinite operators that annihilate \ipo). 
Then any Hamiltonian which is a sum of these operators 
has |V>o) as a ground state. 

The optimum ground state (OGS) approach developed 
in [2l3l4| allows for a systematic construction of such 
ground states. For a given class of Hamiltonians one 
searches for subspaces of the space of interaction param- 
eters where the determination of the ground state can be 
reduced to a local problem. OGS are characterized by 
the fact that they are simultaneously ground states of all 
local interactions contained in the Hamiltonian. It should 
be emphasized that despite this reduction to a local prob- 
lem the properties of the ground state still depend on the 
structure of the lattice. 

Generic examples are simple tensor product states like 
the fully polarized ferromagnet | • • • f). However, the 
OGS approach is more flexible and also allows the con- 
struction of antiferromagnetic states that have a more 
complicated structure. This is possible if the ground 
state of the local Hamiltonian is degenerate. For one- 
dimensional quantum spin chains the generic realization 
of OGS are matrix-product ground state d 5 : 6 : 7 . These 
have been studied systematically for spin-1 0, spin-3/2 
|U and spin-2 p| quantum spin chains revealing the ex- 
istence of various phases, e.g. different types of Haldane 
phases or magnetically ordered states. Here all ground 
state properties, especially correlation functions, can be 
evaluated exactly in closed form. It should be mentioned 
that matrix-product states also play an important role 



for the DMRG-technique****** which can be interpreted 
as variational method based on states of matrix-product 
forir ^ ^^^^^^^^ 

However, the true power of the OGS approach is re- 
vealed by the application to spin models in higher dimen- 
sions d > 2. Technically this generalization is achieved by 
the vertex state model approach**. Here to each vertex 
of a vertex model a quantum state is assigned in addition 
to a Boltzmann weight. This allows for a very efficient 
bookkeeping and provides an elegant way of specifying 
the ground state. Ground state properties like correla- 
tion functions are then determined by the partition func- 
tion of a classical vertex model** on the same lattice. 

This approach has already been applied successfully 
to spin-3/2 systems on the hexagonal lattices* and spin- 
2 systems on the square lattice*!. These models can be 
considered as anisotropic generalizations of the valence- 
bond-solid (VBS) model of Affleck et all** since for a 
generic VBS state the spin S and the coordination num- 
ber z of the lattice are related by S = z/2. 

An important point is the introduction of anisotropic 
interactions. This leads to whole families of solvable 
models and allows to study the parameter-dependence of 
the ground state properties, in contrast to the isotropic 
VBS-type models which are usually solvable only at iso- 
lated points. Indeed in both cases****! a quantum phase 
transition from a disordered to a Neel-ordered phase 
was found. As mentioned above the properties of these 
quantum-critical points are determined by the critical 
behaviour of a classical vertex model on the same lat- 
tice. A similar situation has recently been studied by 
Ardonne et al**V They considered the quantum dimcr 
model introduced by Rokhsar and Kivelson**. For this 
model the ground state can be determined exactly. It 
turns out to be an optimum ground state of a very sim- 
ple form, namely an equal amplitude superposition of 
states connected by the dynamics. On a square lattice 
the Rokhsar-Kivelson model is critical. This allows also 
to study the influence of small perturbations on the crit- 
ical behavious*****. For the Rokhsar-Kivelson model it 
turns out that perturbations can be relevant or irrele- 
vant, depending on the lattice structure. 
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The paper is organized as follows. In Sec. [H] we give 
a rather general introduction to the OGS approach. In 
Sec. lIIII the OGS for the spin-2 case on a hexagonal lattice 
are determined. Here also the properties of these states 
are discussed qualitatively. Sec. IIVI contains concluding 
remarks and an outlook on future work. 



II. OPTIMUM GROUND STATES 

In the following we describe the basic idea of the OGS 
approach which will then be applied to spin-2 models on 
the hexagonal lattice. 

As mentioned before the OGS approach starts from a 
general class of Hamiltonians determined by the physi- 
cal problem and its symmetries. First the most general 
Hamiltonian TL with nearest neighbor interactions con- 
sistent with the following symmetries is constructed: 



1. homogeneity in real space, 



2. parity invariance (i.e. hi 



TABLE I: The symmetry-consistent spin-2 basis states in the 
space of local interactions. Here \m,fh)± denotes the (anti-) 
symmetrization \m, m)± = \m, in) ± \m, m), n = m + fa and 
p = ±1 for symmetric and antisymmetric states. 

/i p name state 
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3. rotational invariance in the xy-plane of spin space, 
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4. spin-flip invariance. 



All of these requirements are quite natural for a two- 
dimensional lattice. The Hamiltonian then has the gen- 
eral structure 



(1) 



{id) 



with nearest neighbor sites 

Table [I] shows the most general basis vectors of the 
space of local interactions for spin S = 2 that are con- 
sistent with the symmetries in terms of Sf-eigenstates 



St\m)i 



m l \m) l 



ith to E {0,±1,±2}. (2) 



of the local spin-2 Hilbert space. It contains seven pa- 
rameters ai , . . . , a-j that correspond to rotations of the 
basis vectors in the respective subspaces that do not vi- 
olate the symmetry requirements. 

The general form of the local interaction hij can then 
be expressed through the projectors on the symmetry- 



consistent basis vectors: 

hij = A 4 (|W4)<1'4| + |f- 4 }(w-4|) + 

A 3 + (K + )(4l + b- 3 )(^3l) + 

A 3 ~(MM + b-3><«-3l) + 
A 2 + 1 (h + 1 )(4ll + l^2l)(^2ll) + 
4l(l«M><«ffll + l«±22>(«±22l) + 

^(]«2 ><«5"l + |w-2>(«r 2 |) + 
Ai(K + i)Rl + l^ 11 )(«iul) + 

Ai(K 2 )(« 1 + 2l + l«il2><^12l) + 

Ari(Ki)(«iil + Kii)(^ul) + 

A^Cte)^! + l u -12)( U -12l) + 

*oi I + A02 k + 2 ) («0 + 2 I + 
Aral^dsX^oal + 

Aoi I %L ) ( V 01 1 + A 02 I %2 > ( V 02 I " 



(3) 



In the following we shall assume for convenience that the 
ground state energy of the local interaction is zero, i.e. 
at least one A is equal to zero. The most general Hamil- 
tonian consistent with the symmetry requirements then 
has 21 free parameters (seven aj and 15 A-parameter, one 
of which is used to fix the local ground state energy), in- 
cluding one trivial scale factor. 
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Now the full Hamiltonian is positive-semidefinitc 
since it is the sum of positive-semidefinite operators . 
Therefore its ground state energy E satisfies E > 0. By 
definition an optimum ground state | -0o ) has Eg = It 
can be characterized by the condition 



hij\ipo) = for all (i, j) 



= 0, (4) 



which means that |i/>o) is simultaneously ground state of 
all local interactions hij . Therefore finite-size corrections 
for the ground state energy vanish exactly for all system 
sizes. 

The simplest examples for OGS are tensor product 
states, e.g. the fully-polarized ferromagnet 



fcrro/ 



n 2 > 



(5) 



The corresponding local ground state is |t>4) = \2, 2) (see 
Table 1). It will become a global ground state if one 
chooses all interaction parameters \ p in J3J as positive 
except for A4 = 0. The parameters 01, . . . , £17 are arbi- 
trary. 

Here we are interested in more complex (antiferromag- 
netic) ground states that are not given by simple tensor 
products. We shall now consider a hexagonal lattice with 
periodic boundary conditions. The specific implementa- 
tion is depicted in Fig. ^ which also shows the numbering 
scheme for the sites. The hexagonal lattice is bipartite 
with sublattices 21, 03. 



' A : ? WW 
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(1,D / 


C.1) 




FIG. 1: Numbering scheme for the sites of the hexagonal 
lattice with sublattices 21 and 23. 




distinguished to which we assign quantum states 



*l M3 \S?(Hi^2,te)), for alii €81, (6) 



! n ) ^, M l 5 f(Mi,M2,M3)) ) for allied, (7) 



i.e. the state vector \S?) is determined by the values \x a 
(a = 1,2,3) of the bond variables belonging to vertex 
i. In the following, fi a is a 2-state variable that can be 
depicted graphically as an arrow pointing in or out of the 
lattice site. 

This general scheme of identifying vertices and quan- 
tum states is rather natural. More specifically in the 
following we use the identification 



Sf = rii + ~ or S? =m-^ 



with 



11 > ■ = 2 l n * (+) ~ n< l ] r 



(8) 



(9) 



where and n\ ' are the number of arrows pointing 
out of and into vertex i, respectively. Therefore all quan- 
tum states generated this way are S^-eigenstates of the 
local spin-2 Hilbert space (see eq. ©)• 

There is still some freedom in the choice of the sign in 
the additive constant eq. JHJ. Since here we are interested 
in ground states that are not macroscopically degenerate 
it will be fixed for each sublattice (see Sec. 1111(1 . For the 
same reason we will not consider bond variables /j, a that 
can take more than two states^ which is possible in prin- 
ciple. The prefactors ce' „ 2 „ 3 are not determined by this 
method and can be chosen in a suitable way restricted 
only by the required symmetries. 

In the next step we define quantum states of the full 
lattice from the single-site states introduced above. As 
in the partition function of a classical vertex model the 
bond variables are summed out. For two sites this implies 



* * * * (10) 

K+K -H 



1. Vertex State Model approach 

On the hexagonal lattice there are two types of sites, 
belonging to the two different sublattices 21, 03 (see 
Fig. Correspondingly two types of vertices can be 



or, explicitly, 



a{a 2 2 \l,l)+a 1 2 al\2,l) = \r' J ) 



1 „ 2 1 



'.1, 2\ 



(11) 



Here the only difference is replacing the product of Boltz- 
mann weights by a tensor product in spin space. This 
procedure is generalized to any number of sites in a 
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straightforward way to construct a quantum state on the 
hexagonal lattice consisting of local spin-2 variables: 



and local quantum states: 



i*>>=£ i® in 

{M \im }l\ I Vie® 




(12) 



where the black vertices have to be taken for i £ 21 and 
the grey ones fori S 23. 

So far this procedure is independent of the Hamiltonian 
under consideration, except for the choice of the local 
Hilbert space. In the final step we now have to specify 
the free parameters in the local interaction © and the 
wavefunction (JSJ, in such a way that the quantum 
state H12fl becomes (i) an eigenstate and (ii) a ground 
state of the Hamiltonian JJJ. This can be achieved by 
imposing the optimum ground state condition 



= hi 




= h 



*3 



= o, 



for any pair of nearest neighbors 



(13) 

and all values of 
fx a for a = 1,2,3,4. Since hij is positive-semidefinite 
this means that the corresponding 2-spin states (|10[) are 
ground states of the local interaction. Therefore, by con- 
struction, the global ground state (|12|l contains only local 
ground states and thus the optimum ground state condi- 
tion n\ip ) = is fulfilled. 



III. GROUND STATES 

In the following we will discuss explicit results for 
ground state wave functions. By specifying the sign of 
additive constant in JHJ we find two different types of 
states. 



A. The Weak Antiferromagnet 

First consider the case where the additive constant in 
iJSJ is fixed to +| for sublattice 21 and to — ^ for sublattice 



®, e.g. 



^ - ai |l} and 



=02 10). The freedom in 



choosing the sign in JSJ can be interpreted as a fourth 2- 
state variable (14, attached to the site (see ©,0) which is 
not summed over in the construction of the global state 
(I12JI . It corresponds to a local spin- 1/2 degree of free- 
dom with magnetization +| on sublattice 21 and — \ on 
sublattice <8. In this sense we can view the construction 
described above as starting from a Neel-type reference 
state with sublattice magnetization | which is modified 
by taking into account the (local) vertex configurations. 
Therefore we have the following identification of vertices 



= a\2) 


b\l) 


b\l) 


^= b\l) 


^= c\Q) 


= c|0) 


c|0) 


*i = b\-i) 


^ = a\-2) 


!^ =ab\-l) 


=ab\-l) 


=ab\-l) 


= c|0) 


jf^= c|0) 


J^= c|0) 


= ab\l) 



They contain three real parameters a, b, c £ K and one 
discrete parameter a = ±1. One of the real parameters 
can always be absorbed into the spectral parameters A 
in © . This choice of parameters is the most general one 
allowed by symmetry requirements. 

Summing out an inner bond variable leads to nine 2- 
spin states which become local ground states in the fol- 
lowing: 

ac|2,0) + fj6 2 |l,l) 
ac|0,-2) + ab 2 \ - 1,-1) 
aab\2,-l) + bc\l,0) 
a6|l,-2) + 6ccr|0,-l) 

bc\l,0) a (14) 
6c| - 1,0) CT 
a 2 |-2,2} + (ife 2 |-l,l} 
ab 2 \ - 1, 1) + c 2 |0, 0) 
ct& 2 |1,-1) + c 2 |0,0) 

As the system should be invariant under time reversal 
the spin-flipped reference state, which corresponds to an 
exchange of the sublattices, should also lead to a ground 
state. This second ground state can be constructed by 
the following identification of vertices and local states: 



= b\i) 


*i = C |o) 


^= c|0) 


*i = C |0) 


6|-1> 


= 6|-1) 


= 6|-1> 


= a |-2) 




^ = c|0) 


J^~ = c\0) 


jf^ = c|0) 


= ab\l) 


^ = cr6|l> 


= cr&ll) 


= o|2) 



As a result the following five additional local ground 
states are obtained: 

ac|0,2) + crfe 2 |l,l) 
ac| - 2,0) +ct6 2 | - 1,-1) 

a 2 |2,-2) + (r6 2 |l,-l) (15) 
crab\ - 1,2} + 6c|0, 1) 
ab\ - 2,1) + bca\ - 1,0). 
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The total number of local ground states is therefore 14. 
The following conditions must be satisfied: 

^2 = ^22 = ^11 = ^12 = ^12 = ^02 = ^03 = Kl = \)2 = 

(16) 

and 

A4, A 3 , A 3 , A 21 ,A n CT , A 01 > 0, (17) 

where a is the discrete parameter appearing in the iden- 
tification of vertices with quantum states. The equations 
(|16fl guarantee that the states are eigenstates whereas the 
inequalities l|17(l are needed to make them ground states. 

The global quantum state l|12f> is the twofold degener- 
ate ground state for the Hamiltonian (|3J satisfying 
above conditions l|l(jfl . (|17fl . The degeneracy is exact for 
all lattice sizes. We do not expect the existence of a 
phase transition within the parameter region where the 
optimum ground state is realized. 

The global ground state is effectively controlled by two 
continuous parameters -, - and one discrete one a ± 1. 
The global Hamiltonian can be tuned by a parameter 
manifold of eight continuous (namely the six parame- 
ters in (|17fl and a/c, b/c) and one discrete parameter a 
with still a free choice of the scale. By construction, the 
ground state is antiferromagnetically ordered. From the 
reference state it inherits a Neel-type order. The sublat- 
tice magnetizations ma = (Sf ) for i G 21 and m<s = (Sf) 
for i £ 03 depend on the parameters -, - and a = ±1 
and satisfy 

m := mg = — m© with < m < 2. (18) 

Therefore in general the sublattices are not fully polar- 
ized (m < 2) and we call the state a weak antiferromag- 
net. For \a\ — > oo one can see from (1141) that the local and 
global ground states are the Neel-states with m = 2. On 
the other hand, for c| — > oo the ground states becomes 
• • ■ 0} and has m — 0. 

B. The Weak Ferromagnet 

It is also possible to start from a different reference 
state by chosing the same sign for both sublattices in 
©. This state is ferromagnetically ordered with mag- 
netization |. Using the same conventions as above, we 
obtain the following identification of vertices with local 
quantum states: 



^ = a\2) 






= b\l) 


c|0) 


= c|0) 


C |o) 


*i = b\-i) 


> y <r = 0-&I-1) 


J^= c|0) 


J^~ = c|0) 


5* = c|0) 


= cr6|l) 


^ = ab\l) 


= ab\l) 


= o|2>. 



These vertices lead to nine local 2-spin states: 

ab\2,l) a 
ac|2,0) +ct6 2 |1,1) 
ac|0,2) + ct& 2 |1,1) 
aab\2,-l) + 6c|l, 0) 
ab\ - l,2) + cr6c|0,l) (19) 

&c[l,0) CT 
c 2 |0,0)+<r6 2 |-l,l) 
c 2 |0,0) + cr6 2 |l,-l) 
bc\ - 1,0) CT , 

where a, b, c are real and a = ±1 is a discrete parameter. 

Again the global ground state is obtained by summing 
over all inner bond variables. The spin-flipped reference 
state with magnetization — h will also lead to a ground 
state if the following additional vertices are used: 





*i = C |o) 


^= c|0) 


c|0> 


*i = b\~i) 


= 6|-1) 


= 6|-1) 


= a \-2) 


^ = a\-2) 


7^ =o-6|-l) 


=<7b\-l) 


^ = o-6|-l) 


= c|0) 


c|0> 


*h ± C |o) 


= o-6|l). 



These vertices generate five additional 2-spin states: 

a6|-2,-l) ff 
ac\ - 2,0) +er6 2 | - 1,-1) 
ac\0, -2) + ab 2 \ - 1,-1) (20) 
aab\ - 2, 1) + 6c| - 1,0) 
a6|l,-2) + er&c|0, -1). 

The 14 states l|19fl and (I20[l are local ground states for 

•\? = ^21 = ^2 = ^11 = ^12 = ^03 = \)2 = ^01 = 

(21) 

and 

A4, A 3 17 , X22, A 12 ° ", Aqi, A02, A 01 > 0. (22) 

Again, the corresponding global states (|12fl are controlled 
by two continuous parameters - , - and one discrete one 
a ± 1. The freedom of chosing the reference state such 
that its magnetization is +^ or — i implies that a priori 
the ground state is exactly twofold degenerate. 

The global state generically has sublattice magnetiza- 
tions — 1 < m<8 < \ < trig < 2 or 1 > m<s > \ > ma > 
—22. It is convenient to introduce the order parameter 
tp = m — i which is the magnetization relative to the 
reference state (with additive constant +1/2). Then we 
have 

3 

^ = V>2t = -V><8 with <?/><-. (23) 
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Note that the magnetization m<s on sublattice 03 can 
indeed become negative so that the global ground state 
exhibits an antiferromagnetic ordering. 

Most interestingly the structure of the ground state 
suggests the occurrence of a quantum phase transition in 
the solvable subspace |J2U, f2*2*)l of the weakly ferromag- 
netic ground state. If the product ab of the parameters 
controlling the local ground states becomes large com- 
pared to the other coefficients in l|19|l , the dominating 2- 
spin states contributing to the global state are ab\2, l) a , 
aab\2,—\) and ab\ — 1,2} (and its spin-flipped counter- 
parts). However, the only global states compatible with 
periodic boundary conditions that can be constructed 
from these local states are product states having |2) on 
sublattice 21 and | — 1) on sublattice 03, or vice versa. This 
implies a fourfold degenerate ground state with sublat- 
tice symmetry in the limit ab — * oc. In this limit the 
order parameter becomes ip = |. On the other hand, 
for small ab the ground state is only twofold degener- 
ate. Thus we can expect a quantum phase transition at 
a critical value (ab) c from a fourfold degenerate state to a 
twofold degenerate state with vanishing order parameter 
■0 = 0, i.e. with non-broken sublattice symmetry. The 
determination of correlation functions requires studying 
a two-dimensional classical vertex model (see Sec. IIVII . 
Preliminary results^ indeed show strong evidence for a 
phase transition line. Numerical results for a representa- 
tive point indicate a continuous transition. 

IV. CONCLUSIONS 

In this paper we have solved a large class of spin- 2 
quantum spin models on the hexagonal lattice. Using 
the optimum ground state concept realized through the 
quantum vertex model approach the explicit form of the 
ground state wave function could be determined. Fo- 
cussing on weakly degenerate ground states two different 
types of states existing in different subspaces have been 
found. In contrast to the generic VBS statesi&2£, coordi- 
nation number z and spin S are not related by S = z/2. 
Some properties of the states that we have constructed 
can be understood by interpreting the spin S — 2 as 



symmetrization of a spin-3/2 and a spin-1/2. For the 
spin S = 3/2 OGS have been constructed previously^. 
In addition we now have a spin-1/2 degree of freedom at 
each site that allows to define different "reference states" . 

The weak antiferromagnet has a twofold degenerate 
ground state characterized by a sublattice magnetization 
ma = — vn<s with < ma < 2. It exists in a subspace 
controlled by eight continuous and one discrete parame- 
ter. The weak ferromagnet is conveniently described by 
the order parameter yb = m — i such that -0a = — ip<B- 
We have argued that it exhibits a quantum phase tran- 
sition from a twofold degenerate disordered phase with 
ip = to a fourfold degenerate phase with broken sub- 
lattice symmetry and ip > 0. 

In future work we will investigate the ground state cor- 
relations which will help to understand the properties of 
the quantum phase transition in more detail. As de- 
scribed ini& expectation values and correlations are de- 
termined by the partition function of a classical vertex 
model on the same lattice. This is a considerable simplifi- 
cation compared to the generic case where e.g. Quantum 
Monte Carlo calculations require the investigation of a 
classical model in dimension d+1. However, in our case 
the classical vertex model has four states per bond and 
is not exactly solvable. Therefore its properties have to 
be studied numerically using classical Monte Carlo sim- 
ulations. 

We want to emphasize that the construction of quan- 
tum states using vertex models is different from the one 
used in so-called quantum vertex models^. There the 
graphical representation of the states in terms of closed 
loops is used to define the full dynamics of a quantum 
model, i.e. the full Hamiltonian. In our case the ver- 
tex model is only used to specify the ground state. The 
identification of vertex and quantum states does not ex- 
tend to excitations. This offers a larger flexibility of the 
approach. 
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